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The  following  report  was  prepared  under  the  AFOSR  RDL  Summer 
Research  Program  under  subcontract  93-102,  F49620-90-C-9076.  It  was 
administered  under  the  direction  of  the  Materials  Directorate,  Wright  Laboratory, 
Air  Force  Materiel  Command,  Wright  Patterson  Air  Force  Base,  Ohio,  with  Dr. 
R.L.  Crane  as  the  Materials  Directorate  Project  Scientist  (WUD).  This  report 
covers  the  doctoral  research  (EKK)  performed  at  the  University  of  Virginia  on 
computer  simulations  of  NLO  biopolymers. 

EKK  is  very  grateful  for  a  doctoral  fellowship  from  the  University  of 
Virginia  Minority  Programs  (funded  by  the  National  Science  Foundation). 


Section  1 


INTRODUCTION 

Photochromic  biopolymers  are  polypeptides  with  substituent  chromophores. 
These  materials  are  being  extensively  investigated  by  the  Air  Force  as  potential 
nonlinear  optical  materials  for  laser  hardening  applications.  Studies  of  the 
structures  of  biologically-synthesizable  materials  that  may  exhibit  NLO  response 
should  contribute  insight  and  understanding  needed  for  the  invention  of  complex 
biomaterials  having  defined  structures  and  controlled  optical  properties.  Helix-to- 
random  coil  transitions  in  polypeptides  are  induced  by  factors  such  as  changes  in 
solution  polarity,  temperature,  and  pH.  Such  transitions  also  occur  when  azo  dyes 
and  spiropyran  chromophores  are  coupled  to  a  polypeptide  backbone  (  1  -  7  ).  The 
polypeptide-chromophore  systems  are  photoresponsive  and  consequently  undergo 
changes  in  their  peptide  backbone  conformation.  This  observation  has  spurred 
research  interest  in  such  photochromic-polypeptide  systems  as  models  for  coupling 
photochemistry  with  backbone  conformational  changes.  They  are  attractive 
systems  because  the  peptide  backbone  can  exist  in  definite  ordered  structures  such 
as  a-helical  or  P-sheet  conformations. 

The  chromophores  investigated  here  are  the  spiropyrans  (  Fig.  1  ).  Current 
applications  of  the  spiropyrans  include  high  resolution  photography,  optical 
devices,  variable  transmission  devices,  photovoltaic  and  holographic  systems  (  8  ). 
In  the  closed  form,  the  indoline  ( left )  and  benzopyran  (  right )  fragments  of  the 
spiropyran  molecule  are  linked  by  a  spiro  carbon  atom  having  tetrahedral  sp^ 
hybridization.  The  two  fragments  are  thus  orthogonal  to  each  other.  Upon 
irradiation  with  light  in  the  UV  region  of  the  electromagnetic  spectrum, 
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spiropyrans  undergo  ring  scission,  converting  from  the  colorless,  closed  form  to 
colored,  open  merocyanine  form  (  Fig.  2  ). 


(a )  closed  form 


(b)  open  form 

Figure  1.  Closed  and  open  forms  of  indoline  spiropyran 


Ring  opening  occurs  by  heterolytic  dissociation  of  the  bond  between  the 
spiro  carbon  and  the  pyran  oxygen  atoms,  followed  by  molecular  rearrangement. 
The  central  carbon  atom  becomes  sp2  hybridized  and  the  indoline  and  benzopyran 
moieties  become  coplanar.  This  facilitates  delocalization  of  the  n  electrons 
between  the  indoline  and  benzopyran  fragments,  giving  rise  to  an  intense 
absorption  band  in  the  visible  portion  of  the  spectrum.  Thermal  energy  or 
absorption  of  visible  light  causes  the  closed,  spiropyran  structure  to  re-form 
slowly. 

Two  variations  of  the  spiropyran  chromophore,  the  indoline  and 
benzothiazoline  derivatives  (Fig.l  and  Fig.  3  respectively)  have  been  investigated. 
Steric  factors  would  be  considerably  different  in  the  benzothiazoline-based 
chromophore  since  it  lacks  the  3 '-methyl  substituents  (  9  )  present  in  the  indoline 
spiropyran.  (Divalent  sulfur  is  substituted  for  the  methylated  in  the  five-member 
ring .) 


Fig.  2.  Reversible  photochromism  in  spiropyran 


When  the  spiropyran  is  appended  to  a  polypeptide,  accumulation  of  charged 
merocyanine  groups  along  the  peptide  chain  can  alter  the  chain  conformation. 
This  in  turn  influences  the  photochromic  behavior  of  the  spiropyran. 
Conformational  calculations  (  7  )  on  spiropyran-modified  poly-L-glutamate  show 
that  the  polypeptide  chain  conformation  is  perturbed  by  the  open  and  closed 
forms  of  the  chromophore,  with  the  open  form  exerting  a  larger  influence  on  the 
peptide  chain. 


Fig.  3.  Benzothiazoline-spiro-benzopyran 


Various  stable  conformations  exist  in  polypeptide  backbones.  These 
include  the  left-handed  and  right-handed  a-helices,  the  extended  |3-sheet,  the  PGI 
and  PGII  helices  found  in  polyglycine,  and  PPI  and  PPII  helices  found  in 
poly  proline.  The  present  molecular  modeling  calculations  focused  on  polypeptide- 
spiropyran  systems  in  which  the  polypeptide  chain  was  initially  in  the  a-helical  or 
^-conformation.  Both  homopolymer  and  copolymer  chains  have  been  considered. 
Poly-L-alanine,  a  homopolymer,  in  both  the  a-helix  and  P-sheet  conformations 
was  studied  in  detail.  Poly-L-alanine  was  chosen  because  it  is  the  simplest 
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polypeptide  found  experimentally  to  adopt  both  the  a-helical  and  |3-sheet 
conformations.  Though  the  most  stable  conformation  of  an  isolated  poly-L- 
alanine  molecule  is  the  a-helix.  X-ray  diffraction  studies  (  10,  11  )  show  that  the 
molecule  can  also  exist  in  an  antiparallel  or  parallel  (3-sheet  conformation  in  fibers 
mechanically  stretched  in  steam  or  treated  with  suitable  solvents. 

The  copolymers  investigated  were  the  hexapeptide  poly(gly-ser-gly-ala-gly- 
ala),  and  the  septipeptide  poly(asp-arg-ieu-ala-ser-tyr-leu).  These  polypeptides  are 
compositionally  representative  of  the  naturally  occurring  amino  acid  sequences  in 
silk  and  wool,  having  respectively,  the  p-sheet  and  a-helical  conformations. 

Previous  work  (  7  )  modeled  the  spiropyran-peptide  system  by  attaching  the 
spiropyran  as  a  pendant  side-group  on  the  polymer.  We  have  investigated  the 
effects  of  alternate  locations  for  the  chromophore  on  the  conformations  of  various 
polypeptides.  For  example,  both  open  and  closed  forms  were  used  as  a  blocking 
group  at  the  amino  terminus  of  the  peptide.  We  also  investigated  systems  in  which 
two  polypeptide  chain  segments  were  attached  to  the  chromophore.  The  structural 
possibilities  for  such  materials  could  be  very  interesting.  The  helical  or  |3-sheet 
segments  could  be  collinear,  and  the  peptide-chromophore  molecules  able  to  order 
easily,  perhaps  into  a  liquid  crystalline  phase.  Furthermore,  the  presence  of  the 
chromophore  could  preclude  crystallization  and  stabilize  the  liquid  crystalline 
phase.  Upon  irradiation  with  light,  formation  of  the  merocyanine  could  cause  the 
a-helix  segments  to  become  noncollinear,  thus  disrupting  the  packing  of  the  liquid 
crystalline  phase.  Not  only  would  the  material  become  colored,  its  opacity  would 
also  change  as  a  result  of  disruption  of  the  molecular  packing.  On  the  other  hand, 
the  packing  of  the  polypeptide  segments  may  substantially  alter  the  ability  of  the 
spiropyran  to  open  and  close,  affecting  its  sensitivity  to  UV  radiation  and  the 
kinetics  of  its  color  change. 


1.1.  Silk  Biopolymers 


Silk-like  proteins  have  elicited  scientific  interest  as  they  possess  inherent 
properties  which  make  them  suitable  for  a  variety  of  applications  including  high 
performance  fibers  (  12  ),  enzyme  immobilization  substrates  (  13  ),  and  polymers 
with  highly  controlled  crystal  morphologies  (  14  ).  Silk  produced  by  the  species 
Bombyx  mod  is  a  fibrous  material  in  which  the  protein  portion  is  composed  of 
fibroin  (  consisting  of  crystalline  fibers  ),  and  a  water-soluble  component  called 
sericin.  Sericin  contains  a  high  percentage  of  polar  amino  acids,  e.g.  serine  and 
aspartic  acid,  whereas  the  fibroin  abounds  in  nonpolar  amino  acids,  e.g.  glycine 
and  alanine.  The  silk  fibroins  have  been  investigated  (  15-27  )  in  detail.  The 
predominant  conformation  observed  is  the  (3-sheet  conformation  (  Figure  4  ),  but 
other  conformations,  for  example,  a-helical,  polyglycine  II,  and  collagen-like 
conformations  are  also  found  (  15  ).  The  p-sheet  conformation  was  also  proposed 
(  18,  1 9  )  for  keratin  in  the  hair  which  had  been  steam-stretched. 

The  antiparallel  (3-structures  which  consist  of  pleated  ^-sheets,  can  be 
divided  into  4  main  groups  on  the  basis  of  their  X-ray  patterns  (  16,  17  ). 

In  group  1,  the  primary  sequence  of  amino  acids  in  Bombyx  mod  is  Gly- 
Ala-Gly-Ala-Gly-[Ser-Gly(Ala-Gly)n]3-Ser-Gly-Ala-Ala-Gly-Tyr,  where  n  is 
usually  2  or  has  an  average  value  of  2.  The  characteristic  X-ray  pattern  shows  an 
axial  period  of  21  A  corresponding  to  the  hexapeptide  sequence  Ser-Gly-Ala-Gly- 
Ala-Gly. 

Group  2  fibroins  are  composed  of  alanine  and  glycine  predominantly  and 
probably  have  a  repeating  sequence  of  the  type  -[(Gly-Ala)x]n-  where  x  =  4,5  or  6. 

Group  3  fibroins  are  composed  of  long  sequences  of  alanyl  residues  and 
feature  an  X-ray  pattern  similar  to  that  of  poIy-L-alanine. 
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Group  4  fibroins  consist  of  alanine  and  glutamic  acid  residues  of  the 
sequence  [Ala-GlnJr,.  Hydrolysis  of  glutamine  in  the  intact  protein  yields  glutamic 
acid. 

In  general,  (3-structures  contain  20-45  %  glycine,  20-50  %  alanine,  and  10- 
40  %  serine.  Thus,  alternating  copolymers  of  these  amino  acids  might  provide 
useful  models  for  studying  the  structure  of  silk.  Poly(L-Ala-Gly)  is  often  used  as 
the  model  for  the  crystalline  fraction  of  Bombyx  mori  fibroin  since  the  two 
crystalline  forms  of  poly(L-Ala-Gly)  are  isomorphous  with  the  two  crystalline 
forms  of  silk. 

In  the  (3-conformation,  the  peptide  chains  form  sheets  with  hydrogen¬ 
bonding  between  adjacent  chains.  The  amino  acid  sidechains  are  accomodated 
between  the  sheets  which  are  held  together  by  van  der  Waals  forces.  The 
structural  features  are  a  chain  repeat  distance  of  ~7  A,  an  interchain  distance  of 
~4.7  A  in  the  direction  of  the  hydrogen  bond,  and  the  N-H  and  C=0  bonds  nearly 
at  right  angles  to  the  chain  axis. 

Figure  5  shows  a  right-handed  a-helical  conformation.  The  a-helical 
conformation  in  synthetic  polypeptides  was  observed  in  the  X-ray  diffraction  of 
poly(7-benzyl-L- glutamate)  (  20,  21  ).  The  characteristic  structural  features  of  the 
a-helix  are  a  pitch  of  5.4  A  and  a  residue  translation  of  1.5  A,  approximating  18 
residues  in  5  turns.  Hydrogen  bonds  are  formed  between  the  i^h  and  (  i-H4  )th 
amino  acid  residues. 

Computational  chemistry  methods  are  essential  tools  for  the  analysis  of  the 
dynamics  and  structure  of  macromolecular  systems,  and  also  the  theoretical 
prediction  of  the  mechanical  properties  of  such  systems.  Conformational  energy 
calculations  on  representative  polypeptide  models  have  been  used  to  study  the 
structural  characteristics  of  collagen  (  22,  23  ),  poly(glycine)  (  24,  25  ), 
poly(alanine)  (  26,  27),  poly(Ala-Gly)  (  27,  28  ),  and  poly(Glu)  (  27  ).  Recently, 
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Pachter  et  al.  (  27  )  used  molecular  simulations  to  predict  the  “spring-like” 
mechanical  response  of  a-helical  poly(Ala)  and  poly(Glu)  that  have  a  reinforcing 
intramolecular  hydrogen  bonding  network.  Also  investigated  was  poly(Gly-Ala)  in 
the  (3-sheet  conformation.  For  the  sake  of  comparison,  calculations  were  extended 
to  include  rigid-rod- like  molecules  containing  an  amide  bond,  i.e.  poly(para- 
phenylene-terephthalamide)  (PPTA).  Results  show  that  the  silk  fibers  exhibit  high 
compressive  strength  and  as  a  result  do  not  kink  in  compression.  In  contrast,  high 
performance  rigid-rod  fibers,  e.g.PPTA,  demonstrate  high  tensile  strength  but  low 
compressive  strength  (  27  ),  hence  do  kink  in  compression.  The  moduli  calculated 
(  27  )  for  extended  poly(Ala)  and  poly(Glu)  chains  were  160  GPa  and  230  GPa 
respectively.  The  response  of  a-helical  poly  (Ala)  chains  to  strain  revealed  no 
kinking  of  the  biopolymer  in  compression.  A  compressive  modulus  of  60  GPa  was 
calculated.  The  results  of  this  study  show  superior  compressive  mechanical 
properties  of  spider  silk-type  polypeptides  over  the  known  high  performance  rigid 
rod  polymers. 

Synthetic  polypeptides,  for  example,  poly(7-benzyl-L-glutamate),  form 
lyotropic  liquid  crystals  in  solution.  The  supramolecular  structure  of  concentrated 
polypeptide  solutions  resembles  that  of  thermotropic,  cholesteric  liquid  crystals 
(20).  .Above  a  certain  concentration,  the  polypeptide  solution  separates  into  a 
birefringent  phase  composed  of  an  anisotropic  arrangement  of  the  polypeptide 
molecules,  and  a  dilute  phase  which  under  an  optical  microscope,  exhibits  a 
spherulitic  structure  typical  of  polymers  crystallized  from  the  bulk  (  21  ). 
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1.2.  Photochromic  Biopolymers 


Photochromic  biopolymers  have  been  extensively  investigated  (1-7,  27, 
28  )  as  materials  with  unique  nonlinear  optical  and  mechanical  properties.  For 
example,  when  thin  films  of  poly(7-benzyl-L-glutamate)  are  doped  with  a 
spiropyran  chromophore,  reversible  changes  occur  in  its  optical  rotatory  power 
(  3  ).  Azo-modified  poly-L-lysine  shows  a  light-induced  conformational  change 
from  an  extended  (3  chain  to  a  random  coil  (  4  ).  A  copolymer  of  |3-benzyl-L- 
aspartate  and  (3-(m-phenylazo)benzyl  L-aspartate  shows  light-induced 
conformational  changes  from  a  left-hand  to  a  right-hand  helix  (  5  ). 
Photomodulation  of  conformation  by  sunlight  in  a  spiropyran-attached  poly-L- 
glutamic  acid  has  been  investigated  (  6  ). 

Figure  6  shows  the  absorption  spectra  of  poly(L-glutamic  acid)  containing 
41  mol  %  of  spiropyran  in  hexafluoro-2-propanol  as  a  function  of  irradiation  and 
dark  adaptation  time.  The  dark  adapted  polymer  shows  a  broad  absorption 
between  550-280  nm,  with  maxima  occurring  at  495,  370  and  312  nm.  Upon 
irradiation  or  exposure  to  sunlight,  the  absorption  in  the  visible  region  of  the 
spectrum  disappears.  Circular  Dichroism  (  CD  )  studies  (  6  )  show  that  the  dark- 
adapted  polymer  exists  as  a  random  coil,  and  the  chromophore  is  in  its  open  form. 
During  light-adaptation,  the  open  form  of  the  chromophore  converts  to  the  closed 
form,  with  a  concurrent  change  in  the  peptide  conformation  from  a  random  coil  to 
an  a-helix.  During  dark-adaptation,  a  slow  helix  to  coil  dark  reaction  occurs  and 
the  chromophore  converts  to  the  colored,  open  form.  A  study  of  the  kinetics  (  2  ) 
of  the  helix-to-coil  dark  reaction  suggests  a  mechanism  which  involves  helix— > 
solvated  helix— >coil  reactions.  Conformational  calculations  (  7  )  on  spiropyran 
modified  poly-L-glutamate  (Fig.  7)  show  that  the  polypeptide  chain  conformation 
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is  altered  by  the  open  and  closed  forms  of  the  chromophore,  with  the  open  form 
exerting  a  larger  influence  on  the  peptide  chain. 

The  stability  of  the  open  form  could  be  influenced  dramatically  by  the 
packing  efficiency  of  the  crystalline  or  liquid-crystalline  state.  Also,  the  site  of 
attachment  of  the  chromophore  could  influence  the  molecular  rearrangements 
which  occur  during  the  open-close  phenomenon. 


Fig.  6.  Absorption  spectra  of  poly(L-gIutaniic  acid)+41  mol%  spiropyran. 
(1,  dark-adapted  solution;  2,  exposure  to  sunlight  or  irradiation).  (6  ) 
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Fig.  7.  Photochromic  behavoir  of  spiropyran-modified  poly-L-glutamate 


Figure  8  shows  two  a-helical  poly-L-alanine  chains  attached  to  the  open 
form  of  the  spiropyran  chromophore  (  designated  as  pA[a]-M-pA[a]  )g^wbile- 
Figure  9  depicts  the  biopolymer  without  the  spiropyran  chromophore  (desig¬ 
nated  as  pA[a]  ).  Figure  10  shows  two  P-poly-L-alanine  chains  attached  to  the 
open  form  of  the  spiropyran  chromophore  (  designated  as  pA[|3]-M-pA[p] ),  while 
Figure  1 1  shows  the  unsbstituted  poly-L-alanine  (  designated  as  pA[P]  ). 
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Fig.  10.  Model  of  pA  [p]-M-pA  [p] 


Fig.  11.  Model  of  pA  [p].  On  the  right  is  shown  the  view  along  the  chain  axis. 


Section  2 


METHODS 

Modeling  calculations  were  made  using  the  SYBYL  Molecular  Modeling 
Software  package  ( version  6.01  )  (  29  ). 

Molecular  mechanics  (  MM  )  was  used  to  generate  and  evaluate 
conformations  of  chromophore-polypeptide  molecules.  Molecular  dynamics  (  MD 
)  was  used  to  evaluate  the  effect  of  a  chromophore  on  the  conformational 
properties  of  the  polypeptide  and  to  address  the  issue  of  molecular  flexibilty  in  the 
system. 

MM  and  MD  calculations  are  based  on  a  classical  description  of  the 
bonding  in  a  molecular  system.  Bond  lengths,  bond  angles,  torsion  angles,  non- 
bonded  interactions,  and  electrostatic  interactions  are  described  in  terms  of  force 
constants  derived  from  a  large  body  of  vibrational  and  crystallographic  data  (  30  ). 
Each  cycle  of  energy  minimization  can  be  thought  of  as  a  step  in  conformational 
space.  The  MM  approach  is  primarily  concerned  with  geometrical  optimization  of 
molecular  conformations  by  finding  the  structural  geometry  that  results  in  the 
lowest  total  potential  energy  (  31,  32  ).  The  total  potential  energy  of  the  molecule 
is  the  summation  of  the  contributions  by  individual  potential  functions. 

=  XEgfj-efch  ^^torsion  ^  ^vdw  ^  ^elec  (  1  ) 

where 

E-stretch  =  bond  stretching  energy  term 
Ebend  =  angle  bending  energy  term 
Etorsion  =  torsional  energy  term 
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Evdw  =  der  Waals  energy  term  (  Leonard  Jones  nonbonded  potential ) 

Egigc  =  electrostatic  energy  term. 

The  relevant  equations  for  each  of  these  energy  terms  are  as  given  below: 

E,,,etch  =  l/2kiKli  -  lo  (  2  ) 

where  kjS  is  the  bond  stretching  force  constant  in  kcal/mol-A^ 

Ij  is  the  length  of  the  ith  bond 

Iq  is  the  equilibrium  length  of  the  ith  bond 

Ebend  =  l/2k,b(e.:  -  e,;0)2  (  3  ) 

where  kjt’  is  the  angle  bending  force  constant  in  kcal/mol-deg2 
9j  is  the  angle  between  two  adjacent  atoms 
is  the  equilibrium  value  for  the  ith  angle 

Etorsion  =  l/2kiKl  +  s/  Isl  (cOStj)  Isl)  )  (  4  ) 

where  kj^  is  the  torsional  barrier  in  kcal/mol 

s  is  the  eclipsed/staggered  rotational  degeneracy  constant 

Evdw  =  kijnb(i/aiji2  -  2/aij6)  (  5  ) 

where  kij"^  is  the  van  der  Waals  constant  in  kcal/mol 

ail  =  r;j/(rj  +  rj),  r;  is  the  van  der  Waals  radius  of  atom  i  (A) 
rjj  =  the  distance  between  atoms  i  and  j 
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^'atoms 

Eeiec  =  332.171  X  Qi  Qj  /  D^ry  (6) 

i=i  j>i 

where  Dy  =  the  value  of  the  dielectric  function  for  atoms  i  and  j 
Qi  =  the  net  atomic  charge  at  atom  i  (  e- ) 

Ty  =  the  distance  between  atoms  i  and  j  ( A ) 

The  SYBYL/BIOPOLYMER  module  was  used  to  model  the  peptide  chain 
conformations  to  be  attached  to  the  indoline  and  benzopyran  components  of  the 
spiropyran  chromophore.  The  biopolymer  modeling  uses  the  concept  of  a 
monomer  and  builds  and  modifies  structures  on  a  residue  by  residue  basis.  The 
polypeptides  were  modeled  as  a-helices  or  |3-sheets  consisting  of  7  amino  acid 
residues.  It  has  been  reported  (  33  )  that  the  minimum  peptide  main-chain  length 
for  a-helix  formation  corresponds  to  7  amino  acid  residues.  The  C-terminus  was 
blocked  by  an  N-methyl  amide  (-HNCH3  )  group,  while  the  chromophore  was 
used  as  a  blocking  group  at  the  N-terminus  of  the  peptide. 

The  conformational  calculations  on  polypeptides  took  into  account 
electrostatic  interactions.  These  were  calculated  by  a  partial  charge  method  (  the 
Del  Re  method  )  which  employs  the  concept  of  localized  bond  orbitals.  The 
Tripos  force  field  treats  hydrogen  bonds  as  non-directional  and  electrostatic  in 
nature.  During  calculations,  hydrogen  bond  energies  are  accounted  for  by  scaling 
the  van  der  Waals  radius  between  the  hydrogen  bond  acceptors  (  N,  O  )  and 
hydrogens  bonded  to  hydrogen  bond  donors  (  N,  O  ).  Structures  were  minimized 
using  MAXIMIN2  which  is  a  full-relaxation  energy  minimization  scheme 
provided  in  SYBYL.  The  default  minimization  procedure  (  Powell  method  )  was 
used.  This  is  an  iterative  method  in  which  the  atomic  coordinates  are  modified 
from  one  iteration  to  the  next  in  order  to  lower  the  energy  of  the  structure.  The 
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method  is  similar  to  the  conjugate  gradient  method,  but  uses  more  advanced  rules 
to  determine  the  descent  direction. 

MD  is  an  extension  of  the  molecular  mechanics  method  of  conformational 
analysis,  employing  the  same  semiempirical  potential  functions  and  associated 
force  field  to  evaluate  atomic  interactions.  MD  calculations  incorporate  kinetic 
(i.e.  thermal  )  energy  to  provide  a  dynamic  view  of  the  system  as  a  function  of 
time  and/or  temperature.  The  atomic  positions  and  velocities  are  iteratively 
computed  in  femtosecond  (  fs  )  time  steps  to  provide  a  description  of  the  natural 
molecular  motions  available  to  the  system  at  the  prescribed  temperature.  In  simple 
terms,  the  MD  calculations  proceed  as  follows.  First  the  force  acting  on  each  atom 
is  computed  (  from  the  derivative  of  the  potential  energy  )  for  all  the  interactions 
(bond  stretches,  bond  angles,  torsions,  van  der  Waals  and  electrostatic)  in  the 
system.  Subsequent  numerical  integration  of  classical  Newtonian  equations  of 
motion  (  F  =  ma,  v  =  at )  using  the  Verlet  (  or  Leapfrog  )  method  (  34  )  yields  the 
individual  atomic  velocities  and  new  atomic  positions.  In  the  Verlet  method,  the 
velocity  is  calculated  at  odd  half  integral  multiples  of  At,  while  the  position  is 
calculated  at  integral  multiples  of  At.  The  velocities  of  the  atoms  are  then  scaled  at 
each  time  step  to  the  specified  temperature  of  the  simulation  using  the  equation: 

SNkgT  =  X  mjVp  (  7  ) 

where  T  is  the  absolute  temperature,  kg  is  Boltzman's  constant,  m^  and  vf  are  the 
mass  and  velocity  of  the  i'*’  atom,  and  N  is  the  total  number  of  atoms.  After  one 
time  step  (  1  fs  )  these  scaled  velocities  bring  the  atoms  to  new  positions  and 
modify  the  forces  they  experience.  The  forces  are  once  again  evaluated  and  the 
procedure  is  repeated  until  the  simulation  accumulates  a  specified  number  of  time 
steps.  A  detailed  description  of  MD  is  given  by  McCammon  (  35  ). 
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MD  was  used  to  simulate  the  natural  molecular  motions  of  the 
chromophore-bound  polypeptides  at  300  K.  Initially  the  system  was  heated  from  0 
K  to  the  target  temperature  of  300  K  in  50  K  increments.  The  total  simulation  time 
employed  was  62,000  fs  including  10,000  fs  each  of  heating  and  equilibration 
stages  of  the  computation.  The  time  step  used  for  each  simulation  was  1  fs. 
Atomic  positions  and  velocities  were  stored  and  data  analysis  was  carried  out  by 
considering  data  at  50  fs  intervals.  Momenta  were  reset  and  nonbonded  interaction 

o 

lists  were  updated  every  25  fs.  The  nonbonded  cutoff  distance  was  8.0  A.  The 
temperature  coupling  factor,  a  parameter  used  to  impose  gradual  temperature 
changes,  was  10.  All  atoms  were  included  explicitly  in  the  simulation.  The 
simulation  was  carried  out  without  solvent,  representing  the  molecule  as  it  might 
exist  and  behave  in  vacuo. 

The  Tripos  force  field  in  SYBYL  version  6.01  was  used  in  the  molecular 
mechanics  and  molecular  dynamics  calculations.  The  force  field  has  been  shown 
(  29  )  to  produce  molecular  geometries  close  to  the  crystal  structures  of  a  variety  of 
organic  molecules  and  biopolymers.  Calculations  using  the  Tripos  force  field  were 
found  (  29  )  to  compare  well  with  the  Amber  force  field,  developed  primarily  for 
biopolymer  calculations. 
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Section  3 

RESULTS  AND  DISCUSSION 

3.1.  Force  Field 

Obtaining  realistic  results  requires  success  on  two  distinct  levels: 
completeness  of  the  model  and  the  accuracy  of  the  force  field  description.  The 
Tripos  Force  field  and  its  associated  potential  energy  functions  were  unmodified. 
To  test  the  reliability  of  the  Tripos  force  field  in  biomolecular  modeling, 
unsubstituted  models  of  poly-L-alanine  in  the  a-helical  and  p-sheet  conformations 
were  built  and  minimized  using  the  Tripos  and  the  Kollman- All-Atom  force  fields. 
The  results  from  each  were  then  compared  with  experimental  values  obtained  from 
X-ray  diffraction  (  18,  19,  21  ).  Tables  1  and  2  show  a  comparison  of  the  results 
from  the  force  field  methods  with  X-ray  diffraction.  The  tables  show  that  the 
Tripos  force  field  is  suitable  for  biopolymer  modeling 

Although  this  force  field  is  parameterized  for  use  with  organic  molecules 
and  biopolymer  molecules,  its  applicabilty  to  the  spiropyran  molecule,  especially 


Table  1  Internal  coordinates  and  chain  repeat  distance  for  an  a-helix 


— 

— 

X-ray 

Kollman 

(p  (deg) 

-58.0 

I - - - 

1  -59.9 

-57.6 

1 

V  (deg) 

-47.0 

i 

I  -44.5 

-47.1 

C  i 

c(A)  i 

27.0 

27.1 
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Table  2.  Internal  coordinates  and  repeat  distance  for  an  isolated  P-sheet 


X-ray 

Kollman 

(t)  (deg) 

-139.0 

-140.8 

-145.8 

139.0 

139.4 

145.0 

c(A) 

7.0 

6.9 

the  merocyanine  is  questionalble.  The  alternating  bond  sequence  (  double-single¬ 
double  bonds  )  (Figure  12  )  linking  the  two  fragments  of  the  merocyanine,  has 
substantial  electron  delocalization.  All  these  binds  have  character  intermediate 
between  single  and  double  bonds,  and  can  potentially  change  with  the 
conformation.  The  Tripos  force  field  describes  these  bonds  as  immutably  single  or 
double;  hence  the  use  of  the  unmodifield  force  field  may  cause  unrealistic  rigidity 
in  this  part  of  the  molecule.  Figure  13  shows  the  torsional  oscillations  of  the 
single  bond  in  merocyanine-substituted  poly-L-alanine.  The  excursions  from  the 
trans  conformation  are  limited  (  ±  10  °C  )  but  still  considerably  larger  (  ±  5  °C  ) 
than  observed  for  the  double  bonds.  Over  longer  simulation  times,  conformational 
transitions  to  gauche  values  would  be  expected.  Since  the  behavior  of  the 
chromophore  may  not  be  representative,  this  analysis  will  focus  on  the  behavior  of 
the  peptide  chains. 


Fig.  12.  Merocyanine,  showing  bonds  linking  the  two  fragments 
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Fig.  13.  Torsional  oscillations  of  the  single  bond  in  merocyanine 

3.2.  Energetics  of  the  Spiropyran-Biopolymer  Models 

Table  3  shows  the  potential  energies  for  the  indoline  spiropyran-poly-L- 
alanine  systems  with  the  peptide  backbone  in  the  helical  or  beta  sheet 
conformation.  Also  included  in  the  table  are  the  energies  of  closed  and  open  forms 
of  spiropyran,  and  unsubstituted  poly-L-alanine  in  a-helical  and  p-sheet 
conformations,  and  having  a  total  length  comparable  to  that  of  each  of  the 
biopolymer-chromophore  systems.  The  large  difference  in  energy  between  the 
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unsubstituted  oc-helical  and  {3-poly-L-alanine  is  due  to  the  fact  that  the  calculations 
were  made  for  an  isolated  chain  of  j3-poly-L-alanine  in  which  the  stabilization 
effect  of  intermolecular  hydrogen  bonding  is  absent.  It  should  be  pointed  out  that 
a  gain  in  stabilization  energy  (  comparable  to  that  of  the  helix  )  of  the  (5- 
conformation  is  normally  achieved  through  intermolecular  hydrogen  bonding 
between  stacked  parallel  or  antiparallel  sheets.  Table  4  shows  a  similar  calculation 
for  the  biopolymer  coupled  to  benzothiazoline-spiro-benzopyran. 

The  sum  of  the  energies  of  the  isolated  molecules  of  poly-L-alanine  and  the 
indoline  spiropyran  is  ~49  kcal/mol,  while  the  energy  of  the  bonded  pA[a]-S-pA 
[a]  or  pA[a]-M-pA[a]  is  65  kcal/mol.  Thus,  upon  derivatizing  a-helical  poly-L- 
alanine  with  the  chromophore,  there  is  a  significant  increase  in  the  energy  of  the 
system  of  about  16  and  32  kcal/mol  for  the  spiropyran-  and  merocyanine- 
substituted  poly-L-alanine  systems  respectively.  The  increase  in  energy  results 
from  some  interference  between  the  helix  and  the  chromophore.  The  energy 
increase  upon  incorporating  the  chromophore  into  the  p-poly-L-alanine  chain  is 
less  than  observed  for  the  helix. 

For  the  benzothiazline  spiropyran-poly-L-alanine  systems,  an  energy 
increase  of  ~30  kcal/mol  results  upon  attaching  the  chromophore  to  a-helical  and 
(3-poly-L-alanine  chains  The  minimum  energy  conformation  is  the  a-helical  poly- 
L-alanine-chromophore  system  irrespective  of  the  form  of  the  chromophore. 

3.3.  Molecular  Dynamics  Results 

Molecular  dynamics  is  used  to  examine  molecular  flexibility  in  the 
biopolymer-chromophore  systems,  and  the  effect  of  mode  of  attachment  of  the 
peptide  chains  on  the  stability  of  the  closed  and  open  forms  of  the  chromophore. 
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Table  3.  PE  ( in  kcal/mol )  of  indoline  spiropyran-poly-L-alanine  models 


Table  4.  PE  (  in  kcal/mol  )  of  benzothiazoline  spiropyran-poly-L-alanine 
models. 
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The  variation  in  the  end-to-end  distance  of  the  chains  attached  to  the 
chromophores  was  studied,  as  well  as  the  angular  variation  in  the  internal 
coordinates  (  (1)  and  \j/  dihedral  angles  )  which  define  the  structure  of  the  peptide. 
The  0  and  \j/  dihedral  angles  measure  the  torsions  about  the  rotationally  permissive 
bonds  in  the  backbone  of  the  amino  acid  residues,  and  hence,  are  the  most 
important  coordinates  in  determining  the  peptide  chain  conformation.  For  a 
peptide  chain  in  the  cx-helical  conformation,  the  values  of  the  0  and  \j/  dihedral 
angles  are  -58°  and  -47°,  respectively,  while  for  the  |3-sheet  the  values  are 
-139°  and  135°.  Molecular  dynamics  simulations  were  carried  out  at  300  K  for  a 
period  of  42,000  fs,  after  10,000  fs  each  of  heating  and  equilibration.  Shown  in 
Figures  14  and  15  are  the  temperature-time  profiles  during  the  MD  simulation  for 
the  a-helical  or  |3-poly-L-alanine-chromophore  systems  investigated. 

3.3.1.  End-to-End  Distance 

The  extent  of  motion  that  occurred  during  the  molecular  dynamics 
simulation  was  observed  by  plotting  the  end-to-end  distance  (  between  the  amide 
nitrogen  of  residue  1  and  the  carbonyl  carbon  of  residue  7  of  the  chain  attached  to 
the  indoline  group  of  the  chromophore  )  as  a  function  of  simulation  time.  Tables  5 
and  6  show  the  statistical  analysis  of  the  end-to-end  distance  of  the  a-helical  and 
P-chains  of  poly-L-alanine  derivatized  with  the  open  and  closed  forms  of  indoline 
spyropyran,  and  the  unsubstituted  poly-L-alanine  in  either  a-  or  P-conformation. 
The  standard  deviation  gives  a  measure  of  the  flexibility  of  the  system,  while  the 
average  distance  measures  the  degree  to  which  the  biopolymer  chain  is  extended  or 
coiled.  The  values  listed  on  the  left  are  for  the  peptide  attached  to  the  indoline 
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Fig.  14.  Temperature-time  profile  for  pA[a]-M-pA[a] 


26 


group,  while  those  on  the  right  are  for  the  benzopyran  group 

Calculations  for  the  unsubstituted  poly-L-alanine  are  made  for  a  model 
composed  of  20  alanine  residues.  The  N-terminus  of  the  chain  was  blocked  by  an 
N-acetyl  group  and  the  C-terminus  by  N-methyl  amide.  The  first  7  (  of  the  20  ) 
alanine  residues  in  the  resulting  model  were  used  to  determine  the  end-to-end 
distance  values  in  the  table.  In  the  chroraophore-biopolymer  systems,  the  N- 
terminus  of  one  chain  was  blocked  by  the  chromophore  and  the  C-terminus  by  an 
N-methyl  amide  group.  This  blocking  could  influence  the  stability  of  the  chain. 
The  stability  of  the  unsubstituted  biopolymer  is  gained  from  the  intramolecular  i^h. 
(i-(-4)'h  hydrogen  bonds  between  the  alanine  residues.  Shown  in  Figures  16  and  17 
are  plots  of  end-to-end  distance  for  the  indoline  spiropyran  attached  to  a-poly- 
L-alanine  and  (3-poly-L-alanine  respectively.  Also  included  in  the  plots  are  the 
data  for  the  unsubstituted  a-helical  and  3-poly-L-alanine. 


Table  5.  Statistical  analysis  of  end-to-end  distance  of  a-poly-L-alanine 


initial  length  of  helix  =  1 1  A 


average  (A) 

12.7,  11.2 

8.1,  11.9 

8.0 

Std  deviation  (A) 

2.7,  2.3 

1.5,  2.0 

1.2 
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Fig.  16.  End-to-end  distance  as  a  function  of  time  (  a-helix  ) 
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Fig.  17.  End-to-end  distance  as  a  function  of  time  ( j3-sheet ) 
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Table  6.  Statistical  analysis  of  end-to-end  distance  of  P-poly-L-alanine 


initial  length  of  p-chain  =  24.0  A 


Results  show  that  there  is  a  significant  deviation  in  the  end-to-end  distance 
of  both  a-  and  13-poly-L-alanine  chains  in  the  presence  of  the  chromophore.  The 
average  lengthof  the  a-helical  chain  attached  to  the  merocyanine  during  the 
simulation  is  slightly  higher  (  ~2  A  )  than  the  starting  value  of  the  helix  and  3-4  A 
longer  than  the  poly-L-alanine  without  a  chromophore.  This  indicates  that  the 
helix  tends  to  uncoil  and  extend  during  the  simulation,  while  the  unsubstituted 
helix  tends  to  coil  up. 

Conversely,  the  extended  j3-chain,  having  an  initial  length  of  ~24  A, 
assumes  a  more  compact  structure  during  the  dynamics.  For  example,  the 
merocyanine-substituted  biopolymer  shows  a  contraction  of  the  chains  during  the 
dynamics,  resulting  in  an  average  distance  of  ~11  A,  a  value  close  to  the  initial 
value  of  an  a-helix.  Substitution  causes  more  collapse  of  the  P-sheet  than  was 
observed  in  poly-L-alanine  by  itself.  Compared  to  the  merocyanine,  the 
spiropyran  forces  the  peptide  chain,  initially  in  either  the  a-helical  or  P-sheet 
conformation,  into  more  compact  structures  as  indicated  by  the  lower  average 
distance. 
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Figure  16  shows  that  the  presence  of  the  merocyanine  causes  the  a-helical 
biopolymer  backbone  to  uncoil  and  extend.  The  spiropyran  does  not  have  much 
effect.  The  number  of  hydrogen  bonds  in  the  a-helical  poly-L-alanine  chain  is 
shown  in  Figure  18.  The  number  of  hydrogen  bonds  decreases  from  the  initial 
(t=0)  value  of  8.  This  allows  the  biopolymer  chain  to  uncoil  and  extend. 

The  MD  simulation  shows  that  a  hydrogen  bond  is  formed  between  the 
oxygen  in  the  NO2  group  on  the  merocyanine  chromophore  and  the  amide 
hydrogen  of  the  first  alanine  residue  in  the  peptide  chain.  This  hydrogen  bond 
remains  intact  during  the  dynamics  simulation  (as  seen  in  Figure  19),  although  the 
length  varies  slightly  from  an  initial  (t=0)  value  of  1.8  A  (  before  the  heating  cycle 
of  the  dynamics  ).  The  hydrogen  bond  stabilizes  the  helix  for  about  12,000  fs  as 
evidenced  in  the  cluster  of  3  points  in  the  a-helical  region  in  Figure  20. 

The  randomization  of  the  |3-sheet  conformation  is  due  to  the  fact  that  MD 
simulations  were  made  for  isolated  j3-chains.  As  shown  in  Figure  17,  the 
spiropyran  chromophore  also  appears  to  enhance  this  randomization,  whereas  the 
merocyanine  has  little  effect.  The  ^-structure  achieves  stabilization  normally  from 
intermolecular  hydrogen  bonding  between  parallel  or  anti-parallel  sheets.  During 
these  MD  simulations  on  the  isolated  chain,  the  |3-sheet  conformation  occasionally 
gains  stability  by  the  formation  of  intramolecular  hydrogen  bonds  between  the  i^h 
and  (i-i-2)th,  or  i‘^  and  (i-t-3)th  amino  acid  residues  as  the  P-sheet  structure 
randomizes.  Figure  21  shows  the  number  of  H-bonds  formed  in  the  P-structures 
during  the  dynamics. 

The  extent  of  helical  backbone  distortion  depends  on  the  attachment  site  on 
the  chromophore.  The  fluctuations  in  the  peptide  chain  attached  to  the  indoline 
and  benzopyran  groups  of  the  chromophore  are  slightly  different.  Table  5  shows 
the  statistical  analysis  of  the  end-to-end  distance  of  the  chains  attached  to  either 
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bonds 


Figure  18.  Number  of  H-bonds  in  a-helical  structures  (both  chains)  during 
simulation  [♦=  pA[a]-M-pA[a],  •  =  pA[a]-S-pA[a],  O  =  pA[a] 
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Figure  20.  Conformation  of  alanine  residues  attached  to  the  indoline  and 
benzopyran  groups  of  merocyanine.  (  +  =  ala  1  at  indoline  side,  •  =  ala  1  at 
benzopyran  side;  numbers  by  the  points  indicate  the  time  step ) 


side  of  the  chromophore.  Standard  deviations  of  2.7  A  and  2.3  A  are  observed  for 
the  chains  attached  to  the  indoline  and  benzopyran  groups  of  merocyanine,  respec¬ 
tively.  Because  the  average  end-to-end  distance  is  closer  to  the  initial  length  of  the 
helix,  the  chain  attached  to  the  benzopyran  group  shows  relatively  more  helical 
character.  A  similar  result  was  found  for  the  benzothiazoline  spiropyran  with 
attached  a-helical  chains.  The  situation  is  reversed  in  the  closed  spiropyran- 
biopolymer  systems.  The  fluctuations  in  the  end-to-end  distance  are  higher  at  the 
benzopyran  ring,  as  shown  in  Table  5. 

Figures  22  and  23  show  the  molecular  conformations  of  the  polypeptide- 
chromophore  systems  at  5,000  fs  intervals  between  30,000  and  60,000  fs.  The 
starting  model  is  shown  at  the  upper  left-hand  comer  while  the  other  seven 
conformations  are  at  30,000  to  60,000  fs  (  in  intervals  of  5000  fs  )  with  time 
increasing  left  to  right  and  top  to  bottom.  The  hydrogen  bonds  are  shown  by 
dashed  lines.  The  sidechains  have  been  omitted  for  the  sake  of  clarity.  The 
starting  and  ending  conformations  of  a-helical  and  |3-poly-L-alanine  spiropyran 
systems  and  unsubstituted  poly-L-alanine  are  shown  in  Figure  24. 

From  the  end-to-end  distance  analysis,  it  can  be  speculated  that  after  the 
chromophore  absorbs  UV  light  to  form  the  merocyanine,  reconversion  to  the  initial 
closed  form  may  be  retarded.  Although  no  experimental  evidence  of  this 
irreversibility  is  available,  the  fact  that  the  merocyanine-bound  poly-L-alanine 
system  asssumes  a  more  helical  character  for  the  attached  polypeptide  chain 
indicates  that  it  is  more  stable  than  the  spiropyran-poly-L-alanine  system.  The 
spiropyran-biopolymer  structure  (  with  reduced  helical  character)  exists  as  a  more 
compact  structure,  hence,  is  expected  to  be  less  stable  than  the  merocyanine 
analog.  Thus,  the  packing  of  the  poly-L-alanine  segments  may  substantially  alter 
the  abiltiy  of  the  spiropyran  to  open  and  close,  affecting  its  sensitivity  to  UV 
radiation  and  the  kinetics  of  the  color  change. 
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3.3.2.  Chromophore-Copolymer  Systems 


Shown  in  Tables  7  and  8  are  the  statistics  of  the  end-to-end  distance  of  two 
types  of  copolymer  chains  attached  to  the  chromophore.  The  values  on  the  left  are 
for  the  chains  attached  to  the  indoline  side,  while  those  on  the  right  are  for  the 
chains  on  the  benzopyran  side.  As  was  done  previously,  the  value  presented  in  the 
table  is  for  the  first  7  of  the  20  residues  in  the  unsubstituted  copolymer.  As  found 
for  unsubstituted  poly-L-alanine,  the  unsubtstituted  copolymer  chain  shows  the 
smallest  deviation,  indicating  that  the  presence  of  the  chromophore  generally 
causes  increased  distortion  of  the  biopolymer  chain.  For  the  a-helical 
septipeptide,  the  merocyanine- substituted  form  shows  an  average  distance  close  to 
the  initial  value  of  the  helix.  As  observed  for  the  homopolymer  chains,  the 
spiropyran  causes  the  copolymer  chains,  initially  in  either  the  a-helical  or  p-sheet 
conformation  to  assume  more  compact  structures,  while  the  merocyanine  leads  to 
larger  end-to-end  distances  than  for  the  unsubstituted  copolymer. 

Table  7.  Statistical  analysis  of  end-to-end  distance  of  a-poly(asp-arg-leu-ala-ser- 
tyr-leu)  (  Co-II[a]  )  chains  attached  to  indoline  chromophore  and  in  the 
unsubstituted  form. 


initial  length  of  helix  =  1 1 .0  A 


Co-II[a]-M-Co-II[a] 

Co-II[a]-S-Co-II[a] 

Co-II[a] 

average  (A) 

12.0,  9.4 

5.6,  11.5 

12.0 

Std  deviation  (A) 

1.6,  1.33 

1.1,  1.0 

0.3 
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Table  8.  Statistical  analysis  of  end-to-end  distance  of  j3-poly(gly-ser-gly-ala-gly- 
ala)  (  Co-l[(3]  )  chains  attached  to  indoline  chromophore  and  in  the  unsubstituted 
form. 

initial  length  of  p-chain  =  20.0  A 


Co-I[13]-M-Co-I[3] 

Co-i[3]-s-Co-i[3] 

average  (A) 

14.0,  11.2 

6.9,  9.20 

10.4 

Std  deviation  (A) 

2.0,  2.6 

1 

1.6,  1.4 

0.6 

The  homopolymer-chromophore  systems  show  a  higher  deviation  in  the 
end-to-end  distance  of  the  chains  than  the  copolymer-chromophore  systems.  As 
observed  for  the  homopolymer  case,  the  extent  of  distortion  of  the  biopolymer 
chain  depends  on  the  site  of  attachment  of  the  chain.  The  merocyanine  results  in 
more  extended  structures  and  the  spiropyran  results  in  more  collapse  for  the 
polypeptide  attached  to  the  indoline  ring  of  the  chromophore.  Comparison  with 
the  homopolymer-chromophore  systems  shows  that  the  effect  of  the  chromophore 
on  the  biopolymer  chains  depends  more  on  the  starting  conformation  of  the 
biopolymer  than  on  the  specific  amino  acid  sequence.  The  additional 
stabilization  effect  from  the  hydrogen  bond  (  formed  between  the  oxygen  in  the 
NO?  on  the  merocyanine  and  the  amide  proton  of  the  first  residue  )  is  minimal  in 
the  copolymer-chromophore  system  as  this  bond,  even  though  formed  initially,  is 
broken  during  the  simulation. 
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3.3.3.  Variation  in  the  (j)  and  Dihedral  Angles 

Shown  in  Tables  9  and  10  are  the  averages  of  the  internal  coordinates,  (j)  and 
w,  of  the  seven  alanine  residues  in  the  a-  and  j3-chains  attached  to  the  indoline  ring 
of  the  chromophore  (  from  30,000  fs  to  31,000  fs  of  simulation  ).  Figures  25  and 
26  show  the  (j),  Xjt  plot  for  each  alanine  residue  starting  in  the  a-helical  and  (3- 
conformations  respectively.  Both  show  large  excursions  of  the  dihedral  angles 
from  their  initial  values,  suggesting  a  randomization  of  these  conformations. 

Table  9.  Average  values  of  the  internal  coordinates,  ([)  and  \j/,  in  an  a-helical  chain 
of  poly-L-alanine-merocyanine  system 


[Numbers  in  parenthesis  show  standard  deviation] 

Table  10.  Average  values  of  the  internal  coordinates,  (j)  and  \j/,  in  a  (3  chain  of 
poly-L-alanine-merocyanine  system 


residue 

# 

1 

2 

3 

4 

5 

6 

7 

b  -139° 

-92(6.0) 

50(11.0) 

76(16.0) 

-104(11.0) 

-79(16.7) 

-105(34.5) 

-116(34.0) 

V/  135° 

162(11.0) 

65(7.0) 

112(11.0) 

46(11.0) 

-50(11.6) 

59(20.3) 

79(26.7) 

[Numbers  in  parenthesis  show  standard  deviation] 
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To  gain  insight  into  the  progression  of  helical  character  along  the  peptide 
backbone,  a  (j),  \|/  plot  for  the  merocyanine-substituted-a-poly-L- alanine  is 
constructed  for  alanine  residues  1,  4  and  7.  Figure  27  shows  the  ij/,  (j)  values  after 
every  5,000  fs.  The  (j),  \j/  values  of  alanine  residue  1  are  close  to  the  ideal  values 
for  an  a-helix.  (  -47°,  -58°  respectively  )  for  about  12,000  fs  of  the  simulation. 
This  results  from  the  stabilization  effect  of  the  hydrogen  bond  which  is  formed 
between  the  oxygen  in  NO2  of  the  merocyanine  and  the  amide  proton  of  the 
alanine  residue. 


Fig.  25.  Average  (f)  and  \(/  dihedral  angles  of  a-helical  alanine  residues 
(  numbers  indicate  residue  number ) 
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3.4.  Molecular  Dynamics  of  Benzothiazoline-spiro-benzopyran-poly-L- 


alanine  Systems 

Tables  11  and  12  show  the  end-to-end  distance  fluctuations  in  the 
biopolymer  chains  attached  to  the  benzothiazoline  spiropyran.  As  observed  for  the 
indoline  spiropyran-biopolymer  systenms,  the  extent  of  chain  distortion  depends 
on  the  attachment  site  on  the  chromophore.  The  P-sheet  conformation  is 


Table  1 1 .  Statistical  analysis  of  end-to-end  distances  of  a-poly-L-alanine  attached 
to  benzothiazoline  spiropyran  (  BZM/BZS  ).  (Values  on  the  left  are  for  the  chains 
attached  to  the  thiazoline  group,  while  those  on  the  right  are  for  the  chains  attached 
to  the  benzopyran  group.) 


o 

initial  length  of  helix  =  1 1  A 


pA[a]-BZM-pA[a] 

pA[a]-BZS-pA[a] 

o 

average  (A) 

11.4,  16.2 

12.1,  12.0 

Std  deviation  (A) 

3.2,  3.1 

2.0,  3.7 

Table  12.  Statistical  analysis  of  end-to-end  distance  of  P-poly-L-alanine  attached 
to  benzothiazoline  spiropyran  (  BZM/BZS  ). 

initial  length  of  P-chain  =  24.0  A 


1 

pA[p]-BZM-pA[p]  1  pA[p]-BZS-pA[p] 

average  (A) 

12.6,  12.0  9.1,  7.0 

Std  deviation  (A) 

4.0,  3.7  2.7,  1.5 
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destabilized  more  than  the  a-helix.  This  result  is  consistent  with  that  obtained  for 
the  peptide-indoline  spiropyran  system.  The  hydrogen  bond  which  is  formed 
between  the  oxygen  in  the  NO2  in  the  merocyanine  form  of  indoline  spiropyran  is 
absent  initially  in  the  benzothiazoline-biopolymer  system.  The  bond  forms  as  the 
.MD  simulation  progresses,  thus  stabilizing  the  biopolymer  chain  to  some  extent. 
The  values  for  the  spiropyran  and  merocyanine  are  more  similar  than  in  the 
indoline  spiropyran-biopolymer  systems.  This  may  be  due  to  the  fact  that  the 
additional  stability  gained  by  the  merocyanine  (  from  the  hydrogen  bond  )  is  less 
than  in  the  indoline-type  systems.  Figure  28  shows  the  end-to-end  distance  plots 
for  a-helical  poly-L-alanine  chain  attached  to  the  thiazoline  ring  of  the 
benzothiazoline  spiropyran  chromophore,  while  Figure  29  shows  a  similar  plot  for 
(3-poly-L-alanine.  The  closed  chromophore-|3-sheet  structure  randomizes  into  a 
more  compact  structure  as  observed  for  the  indoline  spiropyran-biopolymer 
systems.  The  P-sheet  shows  a  tendency  to  assume  a  helical  conformation  by 
forming  i^^-(i-i-2)‘*’  hydrogen  bonds  during  the  simulation,  as  also  observed  for  the 
indoline  spiropyran  systems.  However,  the  number  of  hydrogen  bonds  formed  is 
fewer  than  in  the  indoline  spiropyran  systems  as  shown  in  Figure  30.  Shown  in 
Figures  31  and  32  are  the  molecular  conformations  of  the  poly-L-alanine- 
benzothiazoline  spiropyran  systems  after  every  5,000  fs  simulation.  The  starting 
structures  for  the  a-helical  and  P-poly-L-alanine  systems  are  shown  at  the  top  left. 
The  time  of  the  figures  increases  from  left  to  right  and  top  to  bottom.  The  starting 
and  ending  conformations  of  a-helical  and  p-poly-L-alanine-benzothiazoline- 
spiropyran  systems  are  shown  in  Figure  33. 
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Fig.  29.  End-to-end  distance  as  a  function  of  time 
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Number  of  hydrogen  bonds 


Fig.  30.  Number  of  H-bonds  formed  in  the  P-structures  (both  chains)  during 
simulation.  [♦  =  pA[p]-BZM-pA[p],  •  =  pA[p]-BZS-pA[p],  O  =  pA[p] 
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3.4.1.  Variation  in  ^  and  \|/  Dihedral  Angles 


Figure  34  shows  the  molecular  conformations  of  the  first  alanine  residue  of 
unsubstituted  poly-L-alanine  during  the  dynamics.  Results  from  similar  (t),!}/  plots 
for  the  benzothiazoline  and  indoline  spiropyran-poly-L-alanine  systems  show  that 
both  variations  of  the  spiropyrans  cause  a  significant  deviation  in  the  internal 
coordinates  of  the  biopolymer  chains  (from  the  ideal  (t),\|r  values  of  -58°  and  -47°, 
respectively)  during  the  MD  simulation.  This  implies  that  the  presence  of  the 
chromophore  causes  the  ordered  a-helical  and  {3-sheet  conformations  to  change 
into  a  random  conformation.  Comparison  of  Figure  20  with  the  (t),\(/  plot  in  Figure 
35  for  the  first  alanine  residues  in  the  polypeptide  chain  shows  that  the 
stabilization  effect  of  the  hydrogen  bond  on  the  conformation  is  less  than  observed 
in  the  indoline  spiropyran-biopolymer  system.  The  internal  coordinates  (j),\|i 
remain  close  to  the  ideal  a-helical  values  for  only  5,000  fs,  as  compared  to  12,000 
fs  in  the  indoline  spiropyran-biopolymer  systems.  This  is  due  to  the  fact  that  the 
hydrogen  bond  is  formed  occasionally  during  the  dynamics  as  the  distance 
between  the  oxygen  atom  and  the  amide  proton  gets  closer  to  -1.8  A  (  the  ideal 
value  for  an  0-H  bond  distance  ).  The  diminished  steric  effect  in  the 
benzothiazoline  chromophore  structure,  as  a  result  of  the  absence  of  the  3'-methyl 
substituents,  reduces  the  fluctuations  in  the  chain  attached  to  the  thiazoline  ring  as 
compared  to  the  indoline  ring  (  of  the  indoline  chromophore  ). 
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Figure  35.  Molecular  conformations  of  alanine  residues  attached  to  the 
benzothiazoline  and  benzopyran  groups  of  pA[a]-BZM-pA[a].  (  +  =  ala  1  on 
indoline  side,  •  =  ala  1  on  pyran  side;  numbers  by  the  points  indicate  the  time 
step ) 
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Section  4 


COxNCLUSIONS 

Computational  chemistry  methods  have  been  used  to  explore  the  properties 
of  poly-L-alanine-spiropyran  and  copolymer-spiropyran  systems.  An 
understanding  of  the  molecular  structure  of  these  photochromic  biopolymer 
systems  helps  in  achieving  high  x(2)  or  x(3)  nonlinear  optical  properties. 
Molecular  mechanics  was  used  to  generate  and  evaluate  conformations  for  the 
chromophore-biopolymer  system.  Results  on  the  energetics  of  both  poly-L-alanine 
and  the  copolymer  systems  indicate  that  an  increase  in  energy  of  about  16  to  30 
kcal/mol  occurs  when  the  indoline  or  benzothiazoline  form  of  the  spiropyran 
chromophore  is  coupled  to  6-  or  a-poly-L-alanine.  The  large  difference  in  energy 
between  the  unsubstituted  a-helical  poly-L-alanine  and  poly-L-alanine  in  the  P- 
sheet  conformation  is  due  to  the  fact  that  the  calculations  were  done  on  an  isolated 
chain  of  P-poiy-L-alanine  in  which  the  stabilization  effect  of  intermolecular 
hydrogen  bonding  is  absent.  Molecular  dynamics  simulations  give  insight  into  the 
effects  of  the  chromophore  on  the  conformational  properties  of  the  polypeptide 
and  the  inherent  flexibility  in  the  chromophore-biopolymer  system.  Statistical 
analysis  of  the  end-to-end  distance  during  the  simulations  shows  that  the 
merocyanine  results  in  more  extended  structures  and  the  spiropyran  results  in  more 
collapse  for  the  polypeptide  attached  to  the  indoline  ring  of  the  chromophore.  The 
hydrogen  bond  which  is  formed  between  the  oxygen  (in  the  NO2  group  )  in  the 
merocyanine  chromophore  and  the  amide  proton  of  the  first  alanine  residue  in  the 
peptide  chain  stabilizes  the  a-helix  to  some  extent.  The  P-sheet  structures  achieve 
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stabiltiy  by  forming  intramolecular  hydrogen  bonds  between  the  i'^  and  the  (i+2)*, 
or  i^^’  and  (i+3)*  amino  acid  residues. 

The  effect  of  the  chromophore  on  the  biopolymer  chains  depends  more  on 
the  starting  conformation  of  the  biopolymer  than  on  the  specific  amino  acid 
sequence. 
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Section  5 


FUTURE  WORK 


Future  work  should  be  carried  out  using  semiempirical  molecular  orbital 
calculations  to  explore  the  photochemical  processes,  electron  distributions  as  a 
function  of  substituents  and  conformation,  and  electronic  interactions  between  the 
biopolymer  and  chromophore. 

The  molecular  packing  and  intermolecular  interactions  should  be  analyzed 
in  relation  to  molecular  structure.  Large  layered  assemblies  of  molecules  should 
be  studied  to  examine  end-group  packing  and  to  provide  adequate  boundary 
conditions  for  the  models. 

These  analyses  would  serve  as  a  useful  tool  for  further  prediction  of 
experimental  properties. 
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